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We develop the Hamiltonian theory of axial perturbations around a general time-dependent spher- 
ical background spacetime. Using the fact that the linearized constraints are gauge generators, we 
isolate the physical and unconstrained axial gravitational wave in a Hamiltonian pair of variables. 
Then, switching to a more geometrical description of the system, we construct the only scalar com- 
bination of them. We obtain the well-known Gerlach and Sengupta scalar for axial perturbations, 
with no known equivalent for polar perturbations. The strategy suggested and tested here will be 
applied to the polar case in a separate article. 
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I. INTRODUCTION AND OVERVIEW 



O Perturbation Theory (both hnear and higher order) is one of the most successful tools in General Relativity (GR) . 

It has been used to find the stability properties of a large variety of background solutions like black holes, critical or 
, cosmological solutions. It is also useful to model the evolution of dynamical processes in astrophysical scenarios that 

■ slightly deviate from an exact symmetry, like the oscillations of a static spherical neutron star or a nearly-spherical 
supernova explosion. In particular, it can be used to investigate the emission of gravitational waves in those processes. 

'"^^ . A central problem in GR perturbation theory, inherited from the diffeomorphism invariance of the full theory, is that 
qh' of isolating the physical degrees of freedom from the gauge-dependent information This can be done by imposing 
I , convenient gauge fixing conditions on the perturbations, as Regge and Wheeler [2| originally did in their study of 
^jPj' perturbations of a Schwarzschild black hole. They and later Zerilli Q succeeded in isolating the two physical degrees 
i_i \ of freedom of the gravitational field around spherical vacuum, by taking suitable linear combinations of the remaining 
perturbations and their radial derivatives. These two variables further decouple due to their different properties under 
parity inversion: the Regge- Wheeler variable is axial and the Zerilli variable is polar. 

A more systematic treatment of the gauge freedom in GR perturbation theory was pioneered by Moncrief 4] in his 
Hamiltonian study of the nonspherical perturbations of Schwarzschild. In a Hamiltonian context the four constraints 
— . . obeyed by the twelve dynamical gravitational variables are the generators of the gauge transformations. Moncrief was 
' able to use this information to perform several canonical transformations which reorganized the original six canonical 
^ . pairs of variables into two physical pairs (equivalent to the Regge- Wheeler and Zerilli variables and their canonical 
I ' momenta) and four gauge pairs in which the momenta were constrained to be zero, without any gauge fixing. The 
\ same technique was later applied to other spherical backgrounds with additional symmetries, like Reissner-Nordstrom 

■ 6], Oppenheimer-Snyder [3] or Friedmann-Robertson- Walker Q, but has never been applied to general spherically 
^ ■ symmetric backgrounds, possibly highly time dependent. Hamiltonian perturbation theory has also been recently 

• ^ , revisited in Quantum Gravity with a cosmological background [9[ . A drawback of the Hamiltonian approach is that it 
' is tied to a particular foliation of the background spacetime, and hence the geometric properties of the gauge-invariant 
^ , variables under coordinate transformations involving time are far from obvious. 

■ - - I A Lagrangian formalism was introduced by Gerlach and Sengupta ^lOi] (GS) to study perturbations around generic 
spherical spacetimes, in which the metric perturbation is geometrically split along the decomposition of the 4d manifold 
into the product of a general 2d Lorentzian manifold Af^ with boundary and the unit 2-sphere S*^. This is a 
highly geometrical framework, in which the meaning of the perturbations is transparent, and which also allows the 
construction of gauge-invariant variables. In the axial case it has been possible to isolate the gravitational degree of 
freedom in a single scalar master variable which obeys a wave equation and can be coupled to any kind of matter, both 
in the background and the perturbations. This master scalar and its equation generalize the Regge- Wheeler variable 
and its equation to the axial perturbative problem around spherical symmetry for any reasonable matter model, and 
hence can be considered as the optimal framework for a perturbative study. Unfortunately, in the polar case there 
is not a master scalar valid for a generic spherical background and any matter model, though there are results for 
some particular cases. For instance, a master Zerilli scalar has been introduced by Sarbach and Tiglio for a 
Schwarzschild background, which was later generalized to nonlinear electrodynamics [l^l, around any background 
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solution of the theory. In references [13|, ll4| the gauge-invariant combinations of the stress-energy tensor were also 
included but still on a vacuum background. 

Both approaches to metric perturbation theory are complementary: the Hamiltonian approach offers a better 
framework to handle gauge-invariance, while the Lagrangian approach gives a clearer picture of the geometrical 
structures being perturbed. This Article proposes a combination of both formalisms to construct gauge-invariant 
scalar perturbative variables containing all physically relevant information concerning the gravitational waves. We 
restrict ourselves to spherical backgrounds, but which can be highly dynamical. For definiteness, the dynamics will 
be introduced using a real massless scalar field, but could be done through any other matter model admitting a 
Hamiltonian description. This Article focuses on the axial subset of perturbations, for which the sought solution is 
the Gerlach and Sengupta master scalar [l^l , previously found using the Lagrangian method only. We show how the 
Hamiltonian way allows a more systematic derivation of this object, and how both approaches mutually relate. Most 
important, this Article prepares the path towards a systematic analysis of the polar problem, for which a general 
gauge-invariant master scalar has never been found. Such analysis will be discussed in a second publication. 

The Article is organized as follows. Section presents a brief review of the ADM formalism, establishing the 
notations for 3d and 4d objects, both in the background and first-order perturbations. Section [llll particularizes to a 
spherical background. Section lTVl restricts to axial perturbations and carries out the complete program of scalar gauge- 
invariant construction, as well as establishing the connection between the Hamiltonian and Lagrangian approaches to 
the problem. We conclude in Section |V] with some remarks. 



II. HAMILTONIAN PERTURBATIONS IN GENERAL RELATIVITY 



A. ADM Hamiltonian formalism 



Given the four-dimensional spacetime {M*, '"'g^j,), we introduce a foliation of 3d spacelike slices as level surfaces 
of the time field t{x). The orthogonal vector defines the projected metric ^^^g^^, = '^'ff^j/ + u^u,y on the slices. We 
introduce coordinates (t, x^) adapted to the foliation, and work with three-dimensional objects. Only in the last part 
of this Article we shall use four-dimensional metric variables in order to compare our results with those from the 
Gerlach and Sengupta formalism. Greek and Latin indices denote 4d and 3d tensors respectively. A Icft-superindex 
indicates dimensionality when confusion may arise. 

The 4-metric is decomposed as customary in the lapse function, the shift vector and the 3-metric on the slices 

a-^ ^ -<^'5", /?, ^ g,, ^ <^V,^, (1) 

with inverse 

5'^' = < V + «"'/3'/3^ (2) 

with Latin indices always raised and lowered with g^^ and gij. 

The gravitational dynamical variables in the ADM Hamiltonian formalism are gij and their conjugated momenta: 

ff^ = {g'^K\ - A-^) , Ms = Vdit^, (3) 

where K'^-' is the extrinsic curvature of the foliation hypersurfaces. 

The spacetime will be assumed to contain a dynamical Klein-Gordon field <i>, whose evolution is controlled by the 
action 

Skg = -\j dx^V^9''V'^.^.'^,u (4) 

(5) 



Its canonical momentum has been defined as 



H^-V^'V^*P, (6) 



and denotes the determinant of the 4-metric. The complete action of the system, with coupling constant IGttGjv = 1 
following [4], is given by 



S = Sg + Skg — 



j dt j d\ {n'^g^J,t + n$,t -an- p'n^) . (7) 
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The Lagrange multipliers a and /3* are associated to the constraints 



1 



fJ-g 



(8) 
(9) 



where Dj is the covariant derivative associated to gij. Variation of the action ([7]) with respect to gtj, Hij, and 11 
gives the evolution equations for the corresponding conjugated variables. 



B. Hamiltonian metric perturbations 

Now suppose that the whole system is perturbed at first order. We define the following special notations: 

C = Sa, B' = 5{I3'), (10) 

h,,=5{g,i), f==5{Ii'^), (11) 

= (5$, p = (5n. (12) 

Notice that we perturb the contravariant components of the shift vector because this will give rise to simpler equations, 
even though the comparison with GS variables will be slightly more involved because Pi is better related to the 4-metric 
[seeEq. (P)]. 

Following Taub [lB| and Moncrief [11 we shall obtain the equations for the linear perturbations using the Jacobi 
method of second variations. The idea is that the second variation of the action ([7|), keeping only terms that are 
quadratic on first-order perturbations, gives an action functional for the perturbations, 



(13) 



There are three kinds of terms. First we have kinetic terms, containing time derivatives of hij and ip. Then we have 
the first variations of the constraints that, under a variation of the effective action ([TS]) with respect to and C, 
give the constraints that must be obeyed by the perturbations, 

5{n) = 0, 6{7i^) = 0. (14) 

And finally we have the second variations of the constraints, which are quadratic in the perturbations (ft.^ , p*^, <y3,p), 
and will give the evolution equations for those perturbations. Even though we started with an exact Hamiltonian 
which was a linear combination of constraints, we end up having a quadratic Hamiltonian which does not vanish on 
shell. 

The constraints Ti and Hi in General Relativity are first-class constraints, and hence generators of gauge transfor- 
mations on the constraint surface in phase space. This identifies the gauge orbits, but in general it is not possible 
to separate explicitly the 2 physical degrees of freedom (4 functions) from the 4 gauge variables and the 4 con- 
strained variables in gij and H*-'. The situation in the linearized theory is simpler, but still only highly symmetric 
background scenarios allow the construction of gauge-invariant algebraic combinations of perturbations and their 
derivatives containing the physical information in the linearized approximation. (See [l| for a discussion of the impor- 
tance of symmetry on the algebraic character of the combinations.) One of such background scenarios is a spherically 
symmetric spacetime, as we shall exploit for the rest of this Article. 

For completeness, we provide the expressions for the first variations of the constraints: 

m) = y (n,, - iffyrfi^ [2f^ + 2h\Ii^^ - i/^'^feff-'"^ (15) 



+ -hk + Ms*.**'' + — + - 7;tJ-gh'^'^,t^,j, 

S{n^) = -2Dk{h,jW'' + g,jp''') +W'D,hji+p^,,+Ilip^,. (16) 

See [iBl for intermediate expressions and techniques to compute these expressions. The corresponding expressions for 
their second variations are 

S^{n) = ^{8p^ + 8fil'''G'\h,jhk'' -2h,''hjk) + imjP'' -8P^PJ^ (17) 
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+ ADjhk'^D'W - 4hij{4D'ip + hk^D'<^)D^<^ + lQ{Dih'^ - D'W)Dkhj'' + 4{2Djhik - 3Dkhij)D''h'^]}, 
S^iUi) = 2p^^i + 2p^''iDihjk-2Dkhij)-4hijDkP^''. (18) 



III. SPHERICAL BACKGROUND 



Let us know restrict to a spherically symmetry background M"^ = x S^, where is the unit two-sphere and 
is a two-dimensional Lorentzian manifold with boundary. We shall use arbitrary coordinates = [t, p) on 
and the usual spherical coordinates = {0, 4>) on 5^. Uppercase Latin indices A, B, C, ... denote objects on Af^ and 
lowercase Latin indices a,b,c,... denote objects on S^. The fact that we use arbitrary coordinates on will later 
allow us to keep track of the tensorial character of the different variables. Therefore we do not impose any condition 
on the lapse or shift, apart from being consistent with spherical symmetry. 

The 4-metric can be 2+2 decomposed as 

{ds^)4 = gABix^) dx^dx^ + r^x^) dn^, (19) 

with dfi^ the round metric of the 2-sphere and qab and r being a metric field and a scalar field on M^, respectively. 
We define the vector field va = r~^r^A- 

Using the radial coordinate p explicitly we can write the background spatial 3-metric as 

ids% = a^it, p)dp^ + r^it, p)dQ^^. (20) 

With a spherically symmetric lapse a = a(i, p) and shift vector /?* (/3(i, p), 0, 0) we have 

{ds^)i = {-a^ + a'^f3^)dt^ + 2a^(3dfAp+ {ds^)3 (21) 
= -a^df + a'^{dp + (3dtf +r^dn^, (22) 

which takes the following matricial form, 

The normal vector to the surfaces of constant t is = (—a, 0,0,0) or u'* = a~^(l, — /3, 0,0). Its orthogonal, radial 
vector is n'* = (0, a~^, 0, 0) or = a(/3, 1, 0, 0). In order to work with more geometrical objects we define the following 
frame derivatives that act on any scalar field /: 

a a 

We now derive the background equations, that will be later used to simplify the coefficients of the equations for 
the perturbations. It is convenient to define the following momentum-like variables, which have a definite tensorial 
character with respect to changes of the p coordinate, 

a?W> ^ 2r'^W^ ^ n 

Ui = , n2 = , Hs = — . (25) 

Ms Ms Ms 



We can write the constraints in terms of these spherical variables, 

Ms 



= Hi - n,) - + \ {li^ + $'^) = 0, 



(26) 



a jig r 



- -(r2ni)' + ^n2 + n3$' = o, (27) 
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so that the action is 

■^S = f dt f dp ar^ 2ni^ + 2n2^+n3$t-a— -/3^ . (28) 
The evolution equations can be obtained by simple variation with respect to different variables: 



-[at-ipa)J = ^(Hi-na), (29) 

a 2 

-ir,t-Pr,p) = -^Hi, (30) 
a 2 

^($,t-/3$,p) = n3 (3f) 

i(n,,-,n,,, = M + l-::(^Ji: + i(n>H-*-). (32) 

- (n^.. - /?n2,,) = kiil + Til- n,n,) + ^ - + i (n^ - $'^) , (33) 

a "^2 ar ar 2 \ / 

i,n3,,-,n3,,) = IIiSL±M + (2l!|I. ,34, 

Note that in spherical symmetry the restriction to vacuum, choosing Schwarzschild coordinates (t, r — p), is given by 
n^" = 0, = and $ = 0, that is, fli = 112 = Ha = $ = 0, which simplifies the previous expressions. In particular, 
the constraints and l|?7|) are then trivially obeyed. This is the case studied by Moncrief [J], but here we want to 
analyze the general case. 

IV. AXIAL PERTURBATIONS 
A. Harmonic expansions 

The tensor spherical harmonics form a complete basis on the 2-sphere to expand a tensor field of any rank. Appendix 
\X\ gives the definitions for the harmonics we shall need in this Article, as well as some of their basic properties and 
the relations with the harmonics used by Moncrief Q . See Ref. [IB] for full definitions in the arbitrary rank case and 
further properties. Tensor harmonics can be separated in two groups according to their polarity: there are polar (or 
even, electric or poloidal) harmonics and axial (or odd, magnetic or toroidal) harmonics. In first-order perturbation 
theory around a spherical spacetime polar and axial harmonics decouple and this Article will only deal with the axial 
art of the problem. Following Regge- Wheeler's notation for the metric perturbations and Moncrief 's notations 
for the momentum and shift vector, we expand the perturbative variables in tensor spherical harmonics: 

oo / 

hijdx'dx^ = H H {-'2{hi)Y'dpX';''adx'' + {h2)T XJ^'abdz^dx^), (35) 

l — l m— — l 

oo I 

—p^Jdx'dx^ ^ {-2(Pi)r dp XPa dx" + {p2)T ab dx'^dx''} , (36) 

1=1 m=-l 

Bdx' ^Y.Y. -{hoTXPadx^ (37) 

l — l m— — / 

C = 0, (38) 

p ^ 0, (39) 

if = 0. (40) 

Note that there are no axial perturbations of the 3d scalars a, <& and 11. That means that the scalar field plays no 
role from the perturbative point of view, though the background scalar field is still instrumental to allow for a general 
dynamical spacetime. As we will see, this does not imply any loss of generality. Different {I, m) harmonic components 
also decouple around spherical symmetry, and so from now on we shall drop them from the perturbative variables, 
assuming that we work with a fixed pair of labels at any time. It is important to note that h2, P2 and Hq are scalars 
under changes of the p coordinate, but hi and pi behave as components of a vector. In the language of Id spacetimes. 



6 



the latter are densities of weight +1, to be compensated with metric factors a to convert them into scalars. This will 
become clearer when comparing with the more geometrical GS approach. 

The variables {hi,pi) and (/i2,P2) form two pairs of canonically conjugated variables, whose evolution is partially 
determined by the arbitrary function ho. For example the evolution equations for the variables hi and can be 
easily obtained by perturbation of formula ([3]) after introducing the expansions ([35^40]) 

-(/ii.t-(/?/ji),p) - 2pi+ni/ii + - f^) , (41) 
a a \r'^ J ^ 

-{h2,t- I3h2.p) = 2p2 + (n2-ni)/i2-— . (42) 
a a 

We shall later obtain the evolution equations for more convenient momenta variables. 



B. Gauge-invariant perturbations 

The action functional for the axial perturbations is 

1 (^2^) axial ^ j dtj dx' [p^^ h.,,t - B^^m + ...] (43) 



j dtl^J dp {pihi^t + P2h2,t) + H[ho] 



(44) 



where the dots denote those terms coming from the second variation of the constraints, which we do not need to 
consider in this subsection. The functional H will be defined below in terms of the first variation of the constraint. 
We have also defined 

Pi = Pi, P2 = -tP2, (45) 

a 

where the star stands for complex conjugation and we have defined 

X^\{l-l)l{l + l){l + 2). (46) 



In term of these variables, the perturbed constraint is given by 



^axial _ ^aMg j {r^Pi),P I oP2 ^ ^^211^2 , 2l{l + \) f r^U^hi 
l{l + 1)1 ar^ a r' 



SiHar'^' ^ T^f^ <i + 2^ + X^-^ + "^^^ ^' ( ) !> (47) 



which in turn defines the functional 

H[ho] ^ - J dx^B'Sin,)''"''^ (48) 

'^p{-'-"{^) Pi+^h„p2 + XaIl2^h2-^-^^r'Ui(l^^ /.ij. (49) 

This functional is the generator of gauge transformations, and of course commutes with itself on shell, 

im,HW^j^Jdp[rV^,g-g^J)l^], (50) 

for arbitrary scalar fields / and g. 

Following Moncrief 0] we perform two canonical transformations to separate the gauge-invariant information from 
the pure-gauge content in the canonical pairs (/ii,pi) and (/i2,P2)- The first canonical transformation constructs the 
gauge- invariant combination ki, also a vector component, 

k,=hi + ^('^^ , k2=h2. (51) 
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It induces the following transformation on the momenta: 



TTl =Pi, 

and can be obtained from the generating function 



7r2 = P2 + 



2r2 



F(pi,p2,fci,fe) ^Piki +P2k2 ^-Piy 
In terms of the new variables we can write the first variation of the axial constraint as 



SiHa 



i axial 



Xa^ig I ^TTz _^ ^112^2 ^ 2l{l + 1) 

l{l + 1) I a r2 ar^ 



"2 ^fc2 



(52) 



(53) 



(54) 



which does not contain tti and therefore commutes with ki. That is, ki is gauge- invariant, as we had anticipated. 
This suggests the second canonical transformation: 



with conjugated momenta 



Pi 



TTi -/(? + !) 



Qi = ki, 



Q2 



7r2 



aW2k2 + ^ ^ ^ 



2r' 



fci 



r2 A 

y W 



(55) 

(56) 
(57) 



We can obtain this canonical transformation from the generating function 



F{Pi,P2,ki,k2) = Piki+P2k2+al{l + l) 



ki 



^ --Hi 



fk2 



2a \ r2 



a-i)a + 2) „ ,2 

n2fc2 



(58) 



The first canonical transformation is independent of the dynamical content of the background spacetime, in the sense 
that it does not contain the background momenta Hi, 112, lis. It is actually identical to that of Moncrief 4]. For the 
sake of clarity, we have separated the influence of the dynamical background into the second canonical transformation, 
which trivializes for any static background. 

At this point we have isolated the physical information of the axial metric perturbation in the pair (QijPi) while 
the (Q2, P2) contains the gauge subsystem. P2 is the generator of gauge transformations. 



ma 



, axial 



X. 



2P. 



l{l + l) a 



(59) 



and hence it is gauge-invariant but constrained to vanish. Its conjugated variable Q2 is gauge-dependent and its time 
evolution is determined by the arbitrary function /iq, which can be used to set any desired value for Q2,t- 



C. Evolution equations 

After replacing the new variables and integrating by parts a number of times, we get the following Jacobi action: 

(60) 

where we have defined the first-order quadratic Hamiltonian 

1 



1 ^^2^^ axial ^ / / dp Pl{Ql,t-{/3Qi),p)+P2{Q2,t-l3Q2,p) + 2P2ho-aH^^^ 



H 



(1) = 



ni(FiQi-P2Q2) 



r'^Pi , , r^Hi 



aPl 
'2l{l + l) 



1(1 + 1) 



(/-!)(? + 2) ni(n2-ni; 

2r2 2 



+ ni 



r^P2 



Ql 



(61) 
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The variation of the action with respect to ho gives the constraint that must be obeyed by the perturbations. This 
constraint now takes the simple form 

P2 = 0. (62) 

This constraint is conserved in the evolution since variation with respect to Q2 gives 

{r'P2),t = (/3r2p2),p. (63) 

As P2 is the generator of the gauge transformations, its conjugated variable Q2 is pure gauge. Its evolution equation 
comes from taking the variation of the action respect to P2, 

-{Q2,t-PQ2,p) = -2^-IiiQ2-^(r^Pi+lil + l)^^Qi) . (64) 
a a aA \ ^ / p 

The initial data for Q2 is gauge, and its evolution is fully determined by the free function Hq. In particular it is 
possible to choose Q2 = initially and take ho so that Q2 = at all times. 

We can obtain the physically relevant equations by variation of the action with respect to the variables (Qi,-Pi). 
This gives rise to a system of two coupled second order equations in p-derivatives, whose principal part is, in matricial 
form, 



+ . . . (65) 



the dots denoting lower order terms in p-derivatives of Qi and Pi. We have divided Qi by a to make it a scalar under 
changes of p coordinate. This is a second order in time evolution system, as corresponds to a single wave-like degree 
of freedom, but it apparently has fourth order in p-derivatives for generic values of the background variable Hi. This 
is false because the 2x2 matrix has always vanishing square, and hence the system has third order at most. Actually 
it has second order, as can be checked by taking the matrix to its Jordan canonical form. Defining the combination 

OTT 

L = Pi+l{l + l) — -Qi (66) 
a 

the system (j65p is equivalent to the pair (we now use the dot and prime frame derivatives to simplify the expressions): 

{r^L) = -2A^, (67) 
_ 1 ^V^)Y + gg^(,2x)--(^-^)y + ^)x, (68) 



^2 



2r2 



which can be clearly combined into a single second order equation for L, the sought generalization of the Regge- Wheeler 
equation for dynamical backgrounds. 

When restricting to vacuum the variable rL/\ is the Cunningham-Price-Moncrief master function [7] that obeys 
the Regge- Wheeler equation, though it is not immediately related to the Regge- Wheeler variable. Using the gauge 
r = p, /? = in vacuum we have Hi = and hence rL = rPi, while the Regge- Wheeler variable is Qi/{a?r). We have 
seen that the former is easily generalizable to a dynamical situation as given in (|66p. but not the latter, because it 
would require dividing by Hi, which may vanish. 



D. The master scalar perturbation 

The physical variables Qi and Pi are not scalars in and therefore their values depend upon the foliation we have 
chosen. It is better to describe the gravitational wave using not only gauge-invariant information, but also foliation- 
invariant information, that is, scalars in M^. Studying the geometric properties of those variables under changes 
of foliation is not simple in the 3+1 notation. Following Gerlach and Sengupta [l3|, we change to the M^-adapted 
framework introduced in Section IlIIl in which the foliation is described by the orthonormal frame {u^^n^). This 
allows us to look for scalars on as expressions which are frame independent. The background momenta can be 
rewritten as 

Hi = -2u\a, (69) 
Ha = -2u^VA~2u'^\A, (70) 
= U^^,A- (71) 
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We see that Hi and II3 are essentially time components of vectors in M^. However 112 is a more complicated object. 
In the GS formalism the axial part of the metric perturbation is decomposed in tensor spherical harmonics: 



S (g^^) dx^dx"" = h^^^dx^dx" = ^ {(/^a)?" Xl^bdx^dx'' + Xl^abdx'^dx''} 

l,m 



(72) 



For a given pair (/, m) the vector Ha and scalar h are related in the following way to the original Hamiltonian variables 



hi = ~ihp)GS, 

h2 = 2{h)GS, 

ho = a^{h*)GS- 

The perturbations Ha and h are gauge-dependent, but the following combination is gauge- invariant 

h 



KA = hA-r^[-^ , , 



(73) 
(74) 
(75) 



(76) 



and fully contains the axial information. Therefore it must be given in terms of the gauge-invariant variables Qi, Pi 
and P2: 



= ~Qi, 

a 



K 



Ho 



P2 




r^Pi + 2l{l + 1) 



Those relations can be inverted, giving 

Qi 
Q2 



2ih)GS: 



r'P2 



(/-l)(/ + 2)«*-t-e*^{rV^ [t-'ka]^^} 



(77) 
(78) 

(79) 
(80) 

(81) 
(82) 



We see that the gauge-invariant Qi is the p component of the gauge-invariant vector —ka- Then Q2 is a scalar in 
M^, but it is gauge-dependent. The momentum Pi is the sum of two parts, the first one being a scalar (the curl 
of the vector ka) and the second one being essentially the off-diagonal component of the symmetric tensor W(aKb)- 
Therefore Pi is not a component of a tensor itself. Finally P2 is, apart from a factor aa, the time component of a 
contravariant vector. Again, it is important to stress that these properties are very easy to obtain in GS formalism, 
but not in the original Hamiltonian formalism, where the variables are well adapted to a 3d point of view. 

We want to construct a scalar from a linear combination of Qi and Pi, and perhaps their radial derivatives. We 
already know that Pi is the sum of a scalar plus a non-scalar, so that we have to find whether it is possible to cancel 
that non-scalar part using Qi. It is possible because: 



2{n^u^ + n^u^)vAKB = -2e^^VAKB - 4:{u^VA){n^ kb) 



(83) 



The first term on the r.h.s. is a scalar and the second term is just — 2Hi^. 

Therefore the linear combination L (j66p we defined in the previous section is the scalar we were looking for. It is 
related to the GS scalar variable as 



1 



L 



1(1 + 1) 



H 



GS 



e^^ \r- 



ka\ 



(84) 



This is the most important result of this paper: we have constructed a gauge-invariant variable fully describing the 
physical content of the axial gravitational wave and then we have shown that it is a scalar, so that it is also independent 
of the coordinate system used on the background spacetime. Note that no other independent scalar can be formed as 
a linear combination of Pi and Qi f^^id their derivatives. For the cases when we Hi vanishes, which is equivalent to 
the background gauge condition p = r and /3 = [see Eq. ([50)) ]. the variable Pi is already a scalar. Hence, it is again 
clear that Pi is a more convenient variable than Qi for these cases. 
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Gerlach and Sengupta showed [Tol| that the master variable ([84)1 obeys the following wave equation 



^l^:M±2luas = S^e^^^^iB, (85) 



1^ 



where the bar denotes the covariant derivative on the manifold and "0^ is a gauge-invariant axial perturbation 
of the energy-momentum tensor. This equation, equivalent to the pair (j67H68p for scalar field matter, is valid for all 
background spherical spacetimes and describes the evolution of the axial gravitational wave coupled to any matter 
model. Of course, different matter models will have additional variables and equations coupled to (|85p but we stress 
the fact that both the form of Has and the form of (|85p will remain unchanged. 

The only issue we still need to analyze is whether the vector field Kj^, which could appear in the ipj^ expression or 
in those other equations for the matter variables, can always be reconstructed from the IIgs scalar. An important 
equation in the GS formalism is 

{I - + 2)ka = le^r^VA - eAs(r^nGs)l^, (86) 

whose time component gives ([5^ . after imposing the constraint P2 = 0, and its radial component gives ([S7)) for scalar 
field matter. This equation can be solved algebraically for ka in terms of Has as long as ipA does not contain the 
symmetrized derivative K(^a\b) or higher derivatives of ka- Second and higher derivatives of ka can be ruled out by 
requiring that the matter stress-energy momentum must not contain second derivatives of the metric, because that 
would change the principal part of the Einstein equations. Symmetrized first order derivatives of ka cannot be ruled 
out on physical grounds because perturbation of covariant derivatives of tensor fields may introduce the term 

gaxzal^ya^^^ ^ ^liMEl. (87) 

This is the only possible source of symmetrized derivatives of ka', all other perturbations of Christoffel give either 
Hgs ot undifferentiated ka terms. However, we haven't found any standard matter model in which (j87p appears 
under perturbation of its stress-energy tensor. 

The combination of Hamiltonian gauge methods with the imposition of having a scalar field on has determined 
uniquely the Gerlach and Sengupta scalar Hgs, a variable containing all information on the axial gravitational wave 
and obeying a master wave equation. There are a number of ways of explaining the meaning of this variable (see for 
instance [l7[), but probably the simplest one is given by the expression 

S {RABcd) - -SI±}1 Y eAB ccd r^Hcs, (88) 

where Y is the scalar harmonic. 



V. CONCLUSIONS 



The use of linear perturbation theory avoids the intrinsic nonlinear character of General Relativity, but introduces 
new problems that must always be addressed in some form. A Hamiltonian approach to GR perturbations allows a 
natural discrimination of the gauge from the gauge-invariant information in the problem, but obscures the geometric 
character of the perturbative variables. A complementary Lagrangian approach offers a better geometric understand- 
ing, but gives no clues on how to separate the physical, unconstrained content of the model at hand. This Article 
proposes a combination of those two approaches to arrive at a scalar, gauge-invariant and unconstrained description 
of the linear perturbations in General Relativity. 

Around a spherically symmetric background the axial and polar subsets of perturbations decouple from each other; 
this Article has focused on the axial subset, for which a metric master scalar was already introduced by Gerlach and 
Sengupta ^IQ] using a purely Lagrangian approach. Generalizing Moncrief 's approach for a Schwarzschild background 
we have reobtained this scalar for a general time-dependent background. First, we have isolated the gauge- 
invariant and unconstrained information in a Hamiltonian pair of variables ((3i,Pi), each obeying a first order in 
time evolution equation. Then we have analyzed the geometric character of these two variables, showing that only 
a particular combination of them forms a scalar under transformations on the reduced (under spherical symmetry) 
spacetime. This scalar is, as was expected, the Gerlach and Sengupta master scalar. 

The corresponding master scalar for the polar sector has never been found for a general time-dependent background, 
and there is no known obstruction for its existence. Such a variable would be relevant to study, for example, the 
matching problem through a timelike surface separating two different physical models (like fluid and vacuum at a star 
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surface [18|). Therefore, the open question is whether the same combination of techniques we have used in this article 
can be successfully applied to the polar gravitational wave. This analysis will be presented in a separate publication. 

Many computations in this Article have been performed with the xAct [l9l [20L [2]| framework for Tensor Computer 
Algebra. More precisely, we have used the package xPert for metric perturbation theory around general spacetimes 
to obtain the formulas of Section |IT1 Then, we have expanded these formulas in tensor spherical harmonics with the 
package Harmonics. 
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APPENDIX A: TENSOR SPHERICAL HARMONICS 



It is possible to construct harmonic bases for tensor fields of arbitrary rank on the two-sphere S^. This appendix 
briefly summarizes the construction of those bases for scalars, vectors and symmetric tensors. For the general case 
see [l|. Coordinates on S'^ will be denoted with lowercase Latin letters a,b,c, ... and the round metric will be jab, 
with unit Gaussian curvature. Its associated Lcvi-Civita connection will be denoted with a colon, such that jab-.c — 0. 
Finally, the volume form will be called eab- 

The spherical harmonics YJ™ form a basis for scalar fields on 5^. A vector basis can be constructed from them by 
differentiation, following Regge and Wheeler, as Za = Y-a and Xa = ^abJ^'^Yc, the former being polar and the latter 
axial. Note that the labels {l,m) are always implicitly assumed in the equations. A basis for polar symmetric tensors 
is given by jabY (pure trace) and Zab = Y-ab + '^'^^^ Tab^ (traceless). Axial symmetric tensors can be expanded using 
the basis Xab = {Xa-.b + Xb.,a)/2. 

These tensor harmonics are related to those of Moncrief [4j| as follows: 



[eijijdx^dx-' 


= -2dpXadx°', 


(Al) 


{e2)ijdx^dx'' 


= Xabdx^dx^, 


(A2) 


{fi)ijdx'-dx^ 


= 2dpZadx°', 


(A3) 


(f2)tjdx''dx' 


= dp^Y, 


(A4) 


{h)ijdx^dx^ 


= Yjabdx^dx\ 


(A5) 


{fi)tjdx'-dx^ 


= Ya-.bdx^dx^. 


(A6) 



Under integration over the two-sphere the tensor spherical harmonics form an orthogonal basis, with the following 
normalizations: 

/ dnY*Y' = <5„.5w, (A7) 

dn j''" Zl Z'^ = l{l + l)du'6^m', (A8) 

dnj'''x:xi = l{l + l)du'S,nm', (A9) 

dnrVz:,z;,, - i|±|)|^„,j„„,, (Aio) 

J^jnrV'x^a.x;,, = ll±^su>5,nm'. (All) 
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